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Abstract 

We present a density functional study of the structural and electronic properties of small Cu n 
(n — 1,4) aggregates on defect-free MgO(lOO). The calculations employ a slab geometry with 
periodic boundary conditions, supercells with up to 76 atoms, and include full relaxation of the 
surface layer and of all adsorbed atoms. The preferred adsorption site for a single Cu adatom is 
on top of an oxygen atom. The adsorption energy and Cu-O distance are Es-a = 0.99 eV and 
d S -A = 2.04 A using the Perdew-Wang gradient corrected exchange correlation functional. The 
saddle point for surface diffusion is at the "hollow" site, with a diffusion barrier of around 0.45 eV. 
For the adsorbed copper dimer, two geometries, one parallel and one perpendicular to the surface, 
are very close in energy. For the adsorbed CU3 , a linear configuration is preferred to the triangular 
geometry. As for the tetramer, the most stable adsorbed geometry for CU4 is a rhombus. The 
adsorption energy per Cu atom decreases with increasing the size of the cluster, while the Cu-Cu 
cohesive energy increases, rapidly becoming more important than the adsorption energy. 



I. INTRODUCTION 

The interaction of metallic clusters with supporting metal-oxide surfaces is a subject of great 
current interest, because of the numerous technological applications of these systems, e.g. in the 
field of thin film growth and catalysis Jl]|| . Important objectives of these studies are to understand 
how the atomic and electronic structure of both subsystems are modified through their interaction, 
as well as the properties of the resulting interface. 

In this paper we focus on the adsorption of small Cu n (n — 1,2,3, 4) clusters on the non-polar (100) 
surface of MgO. This surface has been widely investigated both experimentally and theoretically Q. 
At present, all studies agree that the atomic structure of the undefected surface is very close to an 
ideally truncated bulk, even though, from the electronic structure point of view, it is not completely 
clear whether the reduced coordination at the surface gives rise to a reduced ionicity and band gap 
with respect to the bulk. 

Interest in metal clusters has grown remarkably in the past few years. A major question has been 
to understand the dependance of atomic and electronic structures on cluster size and their evolution 
from the small cluster regime to the bulk. In particular, copper clusters have been extensively 
studied, both experimentally and theoretically. As Cu is characterized by a closed d-shell and a 
single valence electron, an important issue has been to investigate similarities and/or differences 
between Cu n and simple alkali-metal clusters. Theoretical work has been done at different levels 
and the properties of small Cu n (n = 1,5) clusters are now quite well understood. 
Studies of Cu clusters and/or overlayers on MgO (100) are not numerous. While most studies 
conclude that the preferred Cu binding site is on top of a surface oxygen and that the Cu-surface 
charge transfer is very small, controversial results have been obtained on the fundamental issue of 
the Cu-surface bond strength. Experimentally, measurements of the initial sticking probability So 
by medium energy ion scattering (MEIS) || have yielded So ~ 0.5 at room temperature, which has 
been interpreted in terms of weak adsorption of Cu adatoms on the MgO surface. More recently, 
however, a substantially larger value of So, ~ 0.82, has been obtained by thermal desorption(TD) 
techniques Theoretically, early Hartree-Fock calculations JhJ found that Cu binds strongly at 
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Mg vacancy sites, while the binding with an undefected surface at the oxygen site is very weak (0.2 
eV). By contrast, more recent Local Density Functional (LDF) calculations found a rather strong 
reactivity of the undefected surface, with an adsorption energy of 1.4 eV for an isolated Cu adatom 
Jilt & t the oxygen site. Similar LDF computations for the Cu/MgO(100) interface Jl2| yielded a 
work of adhesion of 1.0 eV. However, gradient-corrected density functional cluster calculations |l| 
found that the binding of a single Cu atom at the oxygen site is about 0.3 eV. It appears that the 
uncertainty in the value of the adsorption energy is so large that even the character of the Cu-surface 
bonding is unclear. In fact, while adsorption energies of the order of 1 eV or larger suggest a (weak) 
covalent bond, a polarization mechanism should be primarily involved in the case of binding energies 
of the order of 0.2-0.3 eV. 

In this paper we address this issue by means of first principles density functional calculations 
within the framework of the Car-Parrinello approach |l^] . Our calculations employ a slab geometry 
with periodic boundary conditions, and, in order to avoid spurious interactions between periodic 
images of the adatom(s), large supercells (up to 76 atoms) are considered. Moreover, full relaxation 
of the surface layer and of the adsorbed atoms is allowed for (without symmetry constraints). Our 
results show that theoretical predictions of the Cu-MgO binding energy based on density functional 
calculations are quite sensitive to the choice of the exchange-correlation functional. On the basis of 
our most reliable calculations, which use a gradient-corrected functional known as the Generalized 
Gradient Approximation jlf|, the binding energy for a Cu adatom on top of an oxygen of the 
undefected surface is close to 1.0 eV. Moreover analysis of the electronic charge distribution confirms 
the occurrence of a weak covalent bond between the Cu3d/4s states and the 02p orbitals. The barrier 
for the Cu adatom to jump between neighboring binding sites is ~ 0.45 eV, suggesting a substantial 
mobility of the adatom at room temperature. 

Well defined trends show up in the dependence of energetic, structural and electronic properties on 
the size of the adsorbed cluster. In particular we find that the adsorption energy per atom decreases 
with increasing the size of the cluster, while the Cu-Cu cohesive energy increases, and large adsorbed 
clusters are energetically more stable than smaller ones. Thus Cu n clusters do not dissociate, but 
rather tend to maintain their identity when adsorbed on the surface. We also infer that Cu atoms 
which are deposited on the surface at room (or higher) temperature by thin-film growth techniques 
will tend to aggregate and readily form larger clusters. Moreover, calculations for Cus and larger 
adsorbed clusters (to be reported in detail elsewhere jl6| ) indicate a clear preference for 3D over 2D 
geometries. In this context, it may be interesting to remark that experiments Jl7] , p^ | indicate that 
the growth mechanism for Cu on MgO(100) is Stranski-Krastanov, i.e. three-dimensional islands 
are present on top of a full monolayer. 

The layout of this paper is the following. In section [tj] we give the details of our calculations. 
In secti on |lll| , after presenting some convergence tests, we discuss atomic structures and energetics. 
Section |lV| deals with electronic properties. Both these sect ions a re div ided into three main parts. 
We first discuss t he cle an MgO(l OO) surface (sections [II B, and IV A), we then consider the free 
clusters (sections [II C and [IV B), and finally we deal with the adsorption of copper clusters on 
the surface in sections [II D[ and [V C . Both sections [II an d |fv| en d with comments on observed 
trends in the structural/energetical _(|II D 5 ) and electronic (1VC3) properties. A brief summary 
and conclusions are given in section [\J 



II. COMPUTATIONAL DETAILS 

Calculations were carried out within the Car-Parrinello approach using both the Local Density 
(LDA) Jl§,§o| and the Generalized Gradient (GGA) H approximations. For a few calculations, we 
used also the so-called Becke exchange-only functional [ pl[|22| ] , where gradient corrections are applied 
only to the exchange. For the LDA exchange-correlation energy, we employed the Perdew-Zunger 
parametrization |25j] of Ceperley and Alder's electron-gas results. For the GGA exchange- 
correlation energy, we used the functional given in Ref. |l5J. Unless mentioned, the GGA was 
applied perturbatively, i.e. the total energy was calculated with the LDA charge density. This non- 
self-consistent approach to GGA calculations has been shown to yield energies which are usually in 
excellent agreement with fully-self-consistent GGA calculations (see, e.g., Ref. Jl5"|). 

Valence-core interactions were described via Vanderbilt 1251 pseudopotentials for O and Cu and 
Bachelet-Hamann-Schliiter |^6| in the Kleinmann-Bylander ]27[ form for Mg. Pseudopotential cutoff 
radii were 1.5 and 2.0 a.u. for O and Cu respectively. For Cu, valence electrons included the 3d 
and 4s shells. Pseudopotentials were always consistent with the approximation used in the density 
functional calculations, i.e. GGA pseudopotentials for self-consistent GGA and LDA pseudopoten- 
tials otherwise (see e.g. Ref. [^8|). The smooth part of the electronic wave- functions was expanded 
in plane-waves with a cut-off of 16 Ry. Tests using a higher cutoff of 20 Ry were performed (see 
below). A cutoff of 150 Ry was used for the augmented electron density [[HIESj]. As the supercells 
that we used were quite large, the k sampling was restricted to the F point. 



The simulated systems were enclosed in cubic or orthorhombic cells of sizes ranging from 6 to ~ 20 
A with periodic boundary conditions (PBC). In all our calculations, we first performed an electronic 
minimization with steepest-descent and/or damped dynamics algorithms to bring the electrons in 
the ground state corresponding to a given initial atomic configuration. We then relaxed the ions with 
a coupled electronic and ionic damped dynamics. We used the preconditioning scheme of Tassone 
et al. EjJ] to increase the simulation timestep. Timesteps were in the 1.0 — 1.4 ■ 10~ 4 ps range with 
a fictitious electronic mass fi = 1000 a.u.. 

As a first test of the accuracy of our computational scheme, we calculated the equilibrium distance 
(d) and the vibrational frequency (v) of a few molecules relevant to our study, namely Cu2, CuO, 
and MgO. The results, summarized in Table |l| are in very good agreement with the experiment. We 
also determined the equilibrium lattice constant a and the bulk modulus B of bulk MgO. For these 
calculations cubic supercells with 64 atoms were used with only V k-sampling. Also in this case, our 
results for a and B are close to the experimental values (see Table 0) . 



III. STRUCTURE AND ENERGETICS 



A. Convergence tests 

Surfaces were modelled using a repeated slab geometry with periodic boundary conditions parallel 
to the surface. A vacuum of thickness dv is introduced between slabs, dv should be large enough to 
avoid spurious interactions between slabs. Each slab is composed of Nl layers, and surface supercells 
containing (7V at /2) magnesium and (Nat/2) oxygen atoms per layer are used. In our structural 
optimizations, the lower surface of the slab was kept fixed in a bulk-terminated configuration, while 
other layers were fully relaxed. 

Tests with different values of dv, Nl and N a t were performed for a few selected properties. In 
Table III the results of these tests for the surface energy E 3ur f, the binding energy Es-a of a single 
adsorbed copper atom at the on-top oxygen site (see Eq. ^), and the diffusion barrier Ediff of a 
single adsorbed copper atom are reported. E<u //is defined as the difference between the Cu-surface 
binding energies at the absolute minimum (on-top of an oxygen surface atom) and at the saddle 
point (hollow site, see below). E 3ur f is defined as : 

E sur f = — [El°a b — Nl— ^Ebuik] 

ZA ce u Z 

where Elf* b is the total energy of the slab, Ebuik is the energy of an Mg-0 pair in the bulk, and 
A ce u is the area of the supercell. The factor 2 accounts for the two exposed surfaces. Ebuik was 
obtained from a bulk calculation using a supercell with 64 atoms. 

From Table |ll]|, we can see that neither E sur f, Es~a nor Ediff depend significantly on dv ■ 
E sur f depends weakly on Nl, whereas it depends strongly on N a t, i.e. on the size of the supercell. 
This is in turn equivalent to a dependence on the k-sampling of the surface Brillouin zone (larger 
supercells corresponding to improved sampling). The variation of Ediff with Nl and N a t is very 
small. The behavior of Es-a is slightly more complicated, but we can see that also this quantity is 
very stable when large enough supercells are used. 

In Table III results of convergence tests with respect to the plane-wave energy cut-off E™ ut are 
also reported. We note no substantial change bewteen the 16 Ry and the 20 Ry results for Ediff 
and Es-a (in the case N at = 18). Thus, we take E™ u t = 16 Ry as our standard energy cut-off. On 
the basis of these tests, we estimate that our calculated energy differences are accurate within 0.1 
eV. 



B. Clean MgO surface 



Before studying copper adsorption on MgO, we start by a characterization of the clean surface. 
The structure of the fully relaxed surface is very close to the ideal one, as reported in other works 
[ ^o|j3l| . We find an inward surface relaxation of 1.2 % of the bulk inter-layer spacing (0.02 A). This 
is accompanied by a rumpling of 1.5%, with an outward displacement of the oxygens with respect to 
the magnesium atoms. Our best estimate for the surface energy is 1.04 J/m 2 within the LDA (see 
Table |fin|). This value agrees with the results of previous theoretical work | ^o]j32| . The GGA leads 
to a remarkable decrease of the surface energy, resulting in E sur f = 0.86 J/m . A similar effect of 
GGA vs LDA has been found in recent work on Sn02 and Ti02 surfaces |B3|. 



C. Free Cu„ clusters 



Cohesive energies and average distances for free clusters relevant to the present study are reported 
in Table |v| (see also Fig. |l|). The calculations were performed placing the clusters in periodically 
repeated cells of sizes equal to the ones used for the adsorbed systems. The optimized structures of 
these clusters were obtained starting from the geometries of the adsorbed clusters on the surface. 
Thus these structures generally correspond to a local minimum in the potential energy surface. The 
cohesive energy of the free clusters is calculated as: 

D [n] = —Ecu n + nEcu! , 

where Ecu n is the total energy of the cluster with n atoms. The energy of the free copper atom, 
Ecu t , was determined using the same pseudopotential and plane wave cut-off employed for the 
cluster. To account for polarization effects due to the unpaired electron in the atomic 4s state, 
we used the empirical correction AElda-lsd = —0.18 eV(n-\ — nj) 2 [ pif , where n-f (rij) is the 
number of spin up (down) electrons. Our values for Do agree very well with those of other theoretical 
works and experiments (see Table ^). The general trend for Do is to increase as n does. Bonding 
distances are well reproduced. The most stable geometry is the obtuse triangle for n = 3 and the 
rhombus for n — 4 . 



D. Supported clusters 

In this section we present our results for the structure and energetics of supported Cu n clusters 
of sizes n = 1, 2,3 and 4 . The energetics will be characterized in terms of two quantities: 

1. The cluster adsorption energy with respect to the substrate, Es-a ■ 

Es-aH = -[£g£, - E slab - E^]/n 

where BJ 1 ^ is the total energy of the adsorbate-substrate system, E s i a \, is the energy of the clean 
surface, and E(T^ e is the energy of the "free" cluster. As mentioned above, this energy is obtained 
by optimizing the structure of the cluster starting from the geometry of the adsorbed cluster on the 
surface. 

2. The intra-cluster binding energy Ea-a '■ 

E A -a[u] = ~[E a c d u n + (n - l)E slab - n E^ u \]/n 

where Eq^ is the total energy for the surface with one adsorbed adatom. 

The structure of the cluster will be characterized in terms of the average adsorbate-substrate 
distance ds-A , the average copper-copper distance cLa-a and the angle a between an adsorbed 
copper atom, its supporting oxygen atom and the underlying magnesium atom. 

All the results for the following sections are summarized in Table [vj. As the results for a single 
adsorbed copper atom show that the on-top oxygen site is strongly favored (see below), we choose 
to adsorb the cluster atoms at oxygen sites. 

1. Cm 

We considered three possible adsorption sites for a single copper atom on MgO(lOO): on top 
of the oxygen site, on top of the magnesium site and between two oxygen sites (hollow site, see 
Fig. ^|). The preferred site is on top of the oxygen atom (see Table JvJ). In the LDA we find an 
adsorption energy of 1.46 eV and an adsorption distance of 1.89 A . The oxygen atom is attracted 
towards the Cu adatom and moves up by 0.1 A with respect to the other oxygen atoms of the 
topmost layer. The (LDA) interpolated potential energy surface for the Cu adatom is shown in 
Fig. ||. The on-top-0 site is the only minimum, magnesium is a maximum, while the hollow site 
is a saddle point and thus corresponds to the transition state for diffusion. The diffusion barrier 
is Ediff = 0.45 eV. Assuming a simple Arrhenius-like expression for the adatom attempt-to-jump 
frequency F, with a typical prefactor To ~ 10 THz, we can estimate that at room temperature 
F ~ Fo exp(— Ediff /kT) ~ 280 KHz, i.e. a residence time of 3.6 x 10 -6 s. This indicates that 
Cu adatom motion should be readily seen, for instance in Scanning Tunneling Microscopy (STM) 
experiments. Following Ref. , we can approximately define the temperature for which the adatom 
diffusion becomes active (i.e. adatoms jump at least once per second) as 



where Do is the prefactor in the diffusion coefRcent and a (~ 3 A ) is the distance between neighbor- 
ing adsorption sites. When diffusion can be assimilated to a random walk (i.e. for Ediff 3> &b T), 
T is related to D |«|. With T ~ 10 THz, D ~ 2.3 x 10" 3 cm 2 /s, and thus T d ~ 180 K 
(with an error margin of around 20 K, corresponding to an error of one order of magnitude in Do). 
This value of Td has been confirmed to be of the correct order by recent Electron Energy Loss 
Spectroscopy (EELS) experiments 

As the LDA is known to greatly overestimate binding energies, (fully self-consistent) GGA cal- 
culations (also including the optimization of the atomic structure) have been carried out. These 
lead indeed to a reduction of the adsorption energy of about 0.5 eV, so that the resulting value is 
Es-a = 0.99 eV. Within the GGA, also a significant increase of the Cu-O distance with respect to 
the LDA takes place (from 1.89 to 2.01 A ). Note however that GGA calculations which employ the 
LDA charge density and geometry yield a value of the binding energy (0.88 eV) which is very close 
to that given by the full GGA calculations. For consistency with GGA results for larger clusters, 
the latter (non self-consistent) value is reported in Table M. From this Table we can also remark 
that the GGA diffusion barrier is 0.45 eV, i.e. almost equal to that found within the LDA. 

Our LDA adsorption energy, adsorption distance and migration barrier agree well with previous 
LDA calculations [jHJ, whereas there are several important discrepancies between our GGA results 
and the calculations of Pacchioni and Rosch ]l3] |, who used cluster models along with the B-LYP 
[ pl|j38[ exchange-correlation functional. At variance with most available (both experimental and 
theoretical) results, these authors found that the Mg site is slightly more stable than the on-top O 
site, while the hollow site is only ~ 0.1 eV higher than the O site. For the latter they obtained 
Es-a ~ 0.3 eV and ds-A = 2.18 A , i.e. values which are much smaller and much larger than our 
calculated adsorption energy and distance respectively. 

There are two main factors which may contribute to the large discrepancies between our results 
and those of Pacchioni and Rosch. One may be related to the convergence properties of their results 
with respect to cluster size and embedding. The other is the fact that the exchange-correlation 
functionals used in the two calculations are different. 

In an attempt to check the influence of the latter factor, we performed a (fully self-consistent) 
calculation using the Becke exchange-only functional (again including the optimization of the atomic 
structure). We stress that this calculation is only meant to test the dependence of Es-a on the type 
of gradient corrected functional, since it is widely accepted that the Becke-exchange-only functional 
is generally less accurate than the GGA (although for many molecules it yields dissociation energies 
which are in quite good agreement with experiment J3!|). Using the Becke functional, we find 
that the Cu-O adsorption distance becomes 2.09 A (closer to the value of Ref. [^3)), while the 
binding energy is now extremely low, 0.19 eV without spin-polarization correction. This surprising 
result confirms a tendency to underbind of the Becke exchange-only functional (see, e.g., Ref. | p0[ , 
where this tendency is found for hydrogen- bonded systems). It also shows that Cu/MgO is a rather 
"difficult" case, with a strong dependence of the adsorption energy on the choice of the gradient- 
corrected functional. Although no well-defined experimental value of Es-a is available which may 
clearly identify the most appropriate functional for the Cu/MgO system, it is important to remark 
that recent EELS measurements [§7| suggest a value of the Cu diffusion barrier which is very close 
to our result of 0.45 eV, whereas a diffusion barrier of only ~ 0.1 eV is inferred from the calculations 
of Ref. 0. 



2. Cu 2 

We considered three different starting geometries for a copper dimer on MgOQOO): one with the 
two Cu's on nearest neighbor O atoms at distance of about 3 A (A, see Fig. m; the second (B) 
having the copper atoms on second neighbor oxygens at distance 4.2 A , and a magnesium atom in 
between them; the third (S, Fig. ^) with the copper dimer perpendicular to the surface on top of 
an O atom. 

After optimization, the dimer in configuration A is slightly stretched from its free geometry (from 
2.18 A in the free dimer to 2.25 A ) to satisfy bonding with the O atoms, a is 166°. The oxygen 
atoms which support the adatoms move by 0.05 A towards them. The Cu-Cu cohesive energive is 
of order 0.8-0.9 eV (see Table [v|). In configuration B, the two copper atoms cannot move as close 
as they would like to optimize their bond, because of the magnesium atom in between them. The 
Cu-Cu distance is now 2.34 A , and the Cu-Cu cohesive energy is decreased accordingly. The two 
supporting oxygen atoms are displaced from their ideal positions towards the Cu's by rather large 
amounts, ~ 0.14 A . Finally in configuration S, the copper dimer bond length is almost identical 
to that of the free molecule, while the Cu-O distance is the same as that for the adsorbed copper 
adatom. The molecule is slightly tilted in the (110) plane, the O-Cu-Cu angle being 175°. 

Within the LDA, the lowest energy structure is A followed by S (for which the total energy is 0.22 
eV higher than for A) and B (+0.33 eV w.r.t. A). Rather unexpectedly, S is found to be the most 



stable structure within the GGA, even though the difference with A is quite small: O.f 5 eV when the 
GGA is applied non-selfconsistently, and 0.04 eV only when a fully selfconsistent GGA calculation 
(including structural optimization) is performed. This indicates a very delicate balance between 
intradimer cohesion (dominant in the S configuration) and dimer-surface binding (dominant for the 
A structure). 

To investigate this issue further, we sampled the total energy surface as a function of the angle 
(9) between the dimer and the surface, in a plane perpendicular to the surface and passing through 
two neighboring oxygen atoms (the (110) plane). We constrained the value of 9 and allowed all 
the other degrees of freedom to relax. The resulting S-^A energy barrier, occurring at 9 ~ 20°, 
is AS dsmer = 0.13 eV (see Fig. ||). Frequent flips between the S and A structures can be thus 
expected at room temperature. 



3. Cu 3 

The lowest energy structure for the free Cu3 cluster is an obtuse triangle, corresponding to a Jahn- 
Teller distortion of the equilateral triangle ^]-[|. By contrast, we find that the linear arrangement 
on nearest neighbor O sites (configuration A, see Fig. ^) is the preferred adsorption geometry of 
the copper trimer. ft is more bound to the surface and tighter (see the value for ds-A) than the 
triangular arrangement (configuration B). This is due to the presence of a magnesium atom on one 
of the edges of the triangle, which prevents the two copper atoms from binding in an optimal way. 
In the triangular arrangement, all three copper atoms are tilted towards each other: a is 166° for 
the central atom and 162° for the other two. The oxygen atoms move only slightly out of their 
equilibrium positions. In configuration A, the copper atoms at the two ends of the line are tilted 
towards the middle adatom to achieve optimal Cu-Cu distances (a = 159°). Oxygen atoms which 
support copper atoms at the end of the line are attracted to them and move up by 0.08 A while no 
difference can be seen for the middle atom. 



4- Cui 

To determine the most stable arrangement of adsorbed CU4 , two different starting geometries 
were considered: a linear configuration on nearest neighbour O atoms (A), and a square geometry 
with the four Cu's on top of nearest neighbor surface O atoms and a Mg atom at the center of the 
square (B), similar to one considered by Pacchioni and Rosch |l^| (see Fig. ^J). 

We find that upon relaxation the square changes its shape to a rhombus, which, as in the case of 
the free CU4 cluster |]-[|] , becomes the lowest energy structure for the adsorbed tetramer. However, 
the adsorbed rhombus is not perfectly planar : the Cu atoms at the end of the short diagonal 
attract their supporting oxygens by 0.11 A (ds-A = 2.11 A ) with a = 154°, while the other two 
Cu atoms, staying almost on top of their oxygens, lower them with respect to other surface atoms 
by 0.07 A (ds-A = 2.17 A) with a = 175° (see Fig. ^). The adsorption energy per atom of this 
structure is very low (see Table [y|), and close to that found by Pacchioni and Rosch |l3| for CU4 in 
a square geometry. 

Upon relaxation, the linear Cu4 cluster (~ 0.7 eV higher in energy than the rhombus) splits into 
two dimers of length 2.28 A , separated by a distance of 3.28 A (see Fig. [jj). However, these dimers 
are not totally independant, since they show a mirror plane symmetry with respect to the center 
of the line: a — 160° for the extreme atoms and a = 174° for the inner atoms. The oxygen atoms 
supporting the extremity of the line move up by 0.06 A , attracted by the adsorbate atoms, while 
the other two lower their positions by 0.02 A . 



5. Trends 

Competition between adsorption on the surface and cohesion among the atoms in the cluster is 
an important feature of the Cu/MgO system. In the n = 2 case, the dimer (configurations A and B) 
is stretched from its equilibrium distance in the gas-phase to take into account interaction with the 
surface. The balance between intra-dimer and surface-dimer bonding is delicate: in fact full GGA 
calculations yield very close energies for the A and S (standing) geometries. In the n = 3 case, the 
linear structure is preferred to the triangular arrangement, which is more stable in the gas-phase, 
because it allows a better bonding to the substrate. For n = 4, instead, the rhombus, which is by 
far the most stable stable isomer in the gas-phase, is also preferred for the adsorbed cluster, in spite 
of its low adsorption energy. In the n = 4 linear case, the cluster prefers to split into two separate 
dimers in order to optimize its structure and its bonding with the supporting oxygens. 



In Fig. |fj our results for the surface-adsorbate (Es-a, c) and d)) and adsorbate-adsorbate intra- 
cluster (Ea-a, a) and b)) binding energies are summarized. We can remark that the effect of gradient 
corrections is more important on Es-a than on Ea-a, and the absolute difference between LDA and 
GGA results for Es-a tends to reduce with n. With increasing n, the adsorption energy per atom 
Es-a decreases, while Ea-a first decreases (from n — 2 to n = 3) and then increases. Calculations 
for n > 4 clusters JTs| confirm the latter trend. This is an indication that, provided deposition on 
the surface is sufficiently "soft", Cu n clusters will tend not to wet the surface but to retain their 
cluster-like character, particularly in the case of large clusters. Our results also show that larger 
(adsorbed) aggregates are energetically more stable than smaller ones. Thus when Cu atoms are 
deposited on the surface to grow a film, at high enough temperature formation of large aggregates 
by coalescence of small diffusing clusters, particularly monomers, could be observed. 

Calculations for n > 4 jl(| also show that 3D geometries are energetically preferred to 2D ones. 
In particular for n = 5 the lowest energy configuration is a square pyramid (Es-a = 0.43 eV, 
Ea—a ~ 1-08 eV, within GGA). The total energy of this configuration is found to be 0.95 eV lower 
than that of the planar trapezoidal geometry. 



IV. ELECTRONIC PROPERTIES 



A. Clean MgO surface 



The (LDA) density of occupied states (DOS) - obtained by artificially broadening the one-electron 
eigenvalues with gaussians of width 0.15 eV - is displayed in Fig. ^. This shows oxygen 2s (region 
I) and 2p bands (regions II and III) separated into smaller bands. 

Since only occupied states are included, magnesium states should not be seen if one considers 
the MgO (100) surface to be completely ionic and hence the magnesium valence states to be empty. 
Recent calculations have nonetheless shown that some charge stays on the magnesium instead 
of going to the oxygens. This is also the case in our study when one looks at the projected DOS of 
Fig. § 

The overall bandwidth (difference between the top of the 02p and the bottom of the 02s bands) 
is 17.0 eV, while the width of the 02p band is 4.0 eV. Corresponding experimental values by Rossler 
et al. jl2| are 20.0 and 5-6 eV respectively. We also calculated the bandgap between occupied and 
empty states and found E gap = 3.0 eV, against an experimental value of 7.8 eV |43;1 for the bulk. 
The surface band gap has been measured and calculated to be of around 4.5 eV |44[ and 5.0 eV 
p5| respectively. The strong underestimate of the gap is a well-known problem of density functional 
theory. 



B. Free Cu„ clusters 



In order to help understanding the bonding mechanism between the Cu n clusters and the 
MgO(100) surface, in Fig. ^ we show the DOS of the free clusters obtained by optimizing the 
adsorbed structures. Similar DOS for Cu n up to n = 5 have been discussed previously by Jackson 
||. Our monomer's DOS shows the small separation between the 4s state (at ~ -4.5 eV) and the 
3d state (at ~ -5.5 eV). For the dimer, this split-off state on the high energy side of the 3d band is 
missing. This indicates that the dimer is a closed-shell system [^J. This characteristics holds also 
for the CU3 and Cu4 linear clusters, whereas a split-off state on the high-energy side of the DOS 
is present for the triangular CU3 and for the rhombus. The character of the high-energy state is 
mainly s-like, as found by Jackson || and Massobrio et al. || . In analogy with the work of Li et al. 
|fn| , we shall call the split-off state 4s*. 



C. Supported clusters on MgO 



In this section we try to elucidate the bonding between the Cu n clusters and the MgO(OOl) 
surface. To this end we shall examine the DOS for the adsorbed clusters as well as the charge 
density (for the case n — 1). 



1. Cm 



The upper panel of Fig. g shows the (LDA) DOS for a single copper atom at the oxygen and 
magnesium sites of the MgO (100) surface. A split-off state appears on the top of the bands (region 



IV). This state is well inside the bandgap of the clean MgO(lOO) surface, and causes its reduction 
from a value of 3.0 eV for the clean surface to 1.5 eV for the on-top O site. In the magnesium 
site case, the feature corresponding to the 4s* state is smaller. For the O site, prominent states are 
present in region II. These states are responsible for the Cu-surface bonding as shown by Li et al. 
and below. 

In Fig. we present the s,p,d decomposition of the DOS for copper at the stable adsorption 
on-top O site. This decomposition was performed by projecting the plane wave basis onto spherical 
waves centered on the Cu, O, and Mg atoms. Inspection of Fig. H shows that the split-off state 
originates as expected from the Cu4s state. Here no contribution from the d states is present. 
Regions II and III consist mainly 02p states and Cu3d states, the contribution from the Cu4s state 
being small but not negligible especially near the band edges. Region I contains almost exclusively 
02s states. Magnesium states appear in regions II and III, indicating that the ionicity of the surface 
might not be complete. 

To further investigate the character of the bonding, in Fig. [l^ we present plots of the (LDA) 
charge density for the various spectral regions of the DOS. Here, the contours of the integrated 
charge density are shown in a plane perpendicular to the surface and containing only O atoms. For 
technical reasons, the part of the wave functions localized at the cores is missing from the plots. 
In region I, the charge density has the character of s states centered on the oxygens. The charge 
is distorted towards the overlying copper atom. In region II, it clearly indicates the presence of a 
bond through overlap between oxygen 2p and copper 3d/4s states. Region III shows little overlap 
between oxygen and copper orbitals. In region IV the split-off state has anti-bonding character. 
In agreement with the LDA calculations by Li et al. the above analysis of the charge density 
indicates that the bonding between a copper adatom and MgO(OOl) is due to a mixing of the Cu 
3d/4s states and the oxygen 2p band. The bonding part lies in region II, while the antibonding 
counterpart is in region IV. 

The fact that the states in region II are those mostly responsible for the Cu-surface bonding is 
further confirmed by Fig. [ll] which compares the DOS for the n = 1 case calculated using the LDA, 
the GGA, and the Becke exchange-only functional. There is clearly a relationship between Es-a 
and the importance of these bonding states. In fact, the LDA which yields the highest adsorption 
energy also shows the most prominent features in the DOS, while in the "Becke-only" DOS these 
peaks are missing. 



2. CU2,CU S ,CU4 



The DOS for the adsorbed dimer shows few closely spaced states near the top of the 0-2p band 
(see Fig. ||). Furthermore, for the standing dimer (S) the density of states in region II is very similar 
to that of a single copper atom at the on-top-0 site. This indicates that the intra-molecular bond is 
very little perturbed in this case. The closed shell nature of the dimer jl| can explain why it prefers 
to adsorb vertically on the surface (at the GGA level). This can also explain why the adsorption 
energy drops by ~ 0.4 eV from n — 1 to n = 2 (see Fig. |(| c) and d)). At variance with the 
free-cluster case, for n = 3 the linear arrangement (configuration A, see Fig. ^) is preferred to the 
triangular one (configuration B). As mentioned previously, this can be attributed to the fact that for 
configuration B an Mg atom is present in between a pair of Cu atoms, thus hindering their bonding. 
From an electronic point of view, the bonding states in region II are much more prominent for the 
linear case than for the triangular one. We can also remark the presence of two separate 4s* peaks 
for the linear arrangement. For n — 4, the linear tetramer splits into two dimers (configuration A 
of Fig. ^) and displays a DOS similar in shape to that of the dimer. A similarity with the linear 
trimer case can be also noted. For the rhombus a well defined band above the top of the 0-2p band 
is present, while the bonding states in region II are less prominent than for the linear arrangement. 
Indeed the substrate-adsorbate binding energy Es-a for the rhombus is lower than for the linear 
tetramer (see Table ^). Viceversa, Ea-a for the rhombus is substantially larger than for the linear 
case. This is consistent with the fact that for free Cu4 clusters the rhombus is much more stable 
than all other isomers §4|. 



3. Trends 



When n increases, the gap between the 4s* state and the top of the nearest occupied band is 
progressively filled; the formation of an independant band is clearly shown in Fig. ||. Another 
interesting feature is the reduction of the DOS feature corresponding to the bonding states in region 
II. This is a confirmation of the fact that with increasing their size the adsorbed copper clusters 
reduce their ability to bond to the surface. Intracluster bonding becomes dominant and the copper 
agregates thus tend to retain their gas-phase structure. 



V. CONCLUSION 



This paper reports a first principles study of copper microclusters adsorbed on an MgO(fOO) 
surface. At variance with previous calculations where cluster models of rather small size were 
generally used to represent the surface, in our study we use a slab model and periodically repeated 
supercells of large size (up to 36 atoms per layer). Moreover full relaxation, without symmetry 
constraints, of all atoms in the cluster and in the surface layer is allowed for. 

In agreement with previous studies, our results show that the preferred site of adsorption for a 
single Cu adatom on defect-free MgO(lOO) is the oxygen site. The calculated binding energy and 
distance are found to depend significantly on the type of gradient corrections used. According to 
our best estimate, obtained using the GGA, the binding energy is of about 1 eV, suggesting a weak 
covalent bond. This is confirmed by the character of the electronic charge density, which shows the 
mixing of copper 3d and 4s states with oxygen 2p orbitals. The Cu hopping barrier is ~ 0.45 eV 
(largely independent of the functional), suggesting that significant adatom diffusion should be seen 
starting from ~ f80 K. 

When increasing the size of the clusters, we find that the intra-cluster binding energy tends to 
dominate over the adsorption energy, larger (adsorbed) aggregates being energetically more stable 
than smaller ones. These results suggest that the clusters which are "softly" deposited on the 
surface at low temperatures will tend to keep their gas-phase structure (see, e.g., Cu2 and CU4) 
while optimizing their bonding with the substrate. Viceversa, when Cu atoms are deposited on 
the surface to grow a Cu film on MgO(fOO), formation of large aggregates by coalescence of small 
diffusing clusters (particulary monomers) should be observed at high enough temperature. Results 
for n > 4 (to be presented elsewhere) show that bigger copper clusters tend to form three-dimensional 
structures instead of planar ones. 
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TABLE I. Calculated (LDA) bond distances and stretching frequencies of a few selected molecules. 



Molecule 


LDA 


d[A] 




Exp. @| 




LDA 


v [eV] 

Exp. 0,0 


Cu 2 
CuO 
MgO 


2.18 
1.69 
1.77 






2.22 
1.72 
1.75 




0.034 
0.087 
0.103 


0.033 
0.079 
0.097 


TABLE II. Calculated structural parameters for bulk MgO. a 
parameters from the Murnhagen equation of state. 


is the lattice constant, B is the bulk modulus, Vo and B' are 












a [A] 


B [Mbar] 


Vo [A 3 ] B' 


This work 


(LDA) 
(LDA) 


Murnhagen fit 
Polynomial fit 


4.25 
4.25 


1.63 
1.57 


68.74 3.70 
68.63 


Theory 


(Hartree-Fock) 


Ref. 


18 




4.20 


1.86 


66.18 3.53 


Exp. 




Ref. 




-|52] 


4.21 


1.55-1.62 


66.60 



TABLE III. Convergence tests for the surface energy (E sur f), the Cu-surface binding energy (Es-a) and the Cu diffusion 
barrier (Ediff). Calculations are within the LDA. Nl is the number of layers, N a t is the number of atoms per layer, d v is the 
width of the vacuum. E™ ut is the energy cut-off for the smooth part of the wave-functions. 



N L 


N at 


Ec Ut [Ry] 


d„ [A] 


Esurf [J/m 2 ] 


Es-a [eV] 


Ediff [eV] 


2 


8 


16 


7.4 


1.59 


1.71 




3 


8 


16 


7.8 


1.60 


1.29 


0.53 


3 


8 


16 


19.8 


1.77 


1.25 


0.51 


3 


8 


20 


13.8 




1.50 


0.53 


2 


18 


16 


12.8 


1.19 


1.46 


0.45 


2 


18 


20 


12.8 




1.50 




3 


18 


16 


10.7 


1.19 


1.48 




2 


36 


16 


12.8 


1.04 







TABLE IV. Cohesive energy per atom (Do/n) and average Cu-Cu distance for the free clusters shown in Fig. [l|. The 
values for the cohesive energy have been obtained within the CCA (LDA). Spin-polarization effects are included in the 
works of Jackson and Calaminici et al. , but not in that of Lammers et al. . 



11 


Geometry 


Ref. 


Do/n [eV] 


d[A] 


2 


Line 


This work 


1.13 (1.33) 


2.18 


2 


Line 


Jackson 


1.08 (1.36) 


2.18 


2 


Line 


Calaminici et al. 


1.13 (1.30) 


2.20 


2 


Line 


Exp 0g 


1.04 


2.22 


3 


Obtuse triangle (dv) 


This work 


1.13 (1.43) 


2.31 


3 


Obtuse triangle (dv) 


Jackson f| 


1.16 (1.52) 


2.27 


3 


Obtuse triangle (dv) 


Calaminici et al. 


1.12 (1.34) 


2.34 


3 


Obtuse triangle (dv) 


Exp. @ 


1.02 




3 


Line 


This work 


1.10 (1.36) 


2.23 


3 


Line 


Lammers et al. M, 


0.57 




3 


Line 


Calaminici et al. H 


1.18 (1.48) 


2.27 


4 


Rhombus (Z?2h) 


This work 


1.53 (1.89) 


2.32 


4 


Rhombus (D 2 h) 


Jackson 


1.52 (1.95) 


2.30 


4 


Rhombus (D^h) 


Lammers et al. H] 


1.04 




4 


Rhombus (D^u) 


Calaminici et al. [H 


1.59 (1.90) 


2.36 


4 


Line 


This work 


1.48 (1.74) 


2.24 



TABLE V. Calculated energies (in eV) and structures (distances in A) for adsorbed copper clusters shown in Fig. ^[ All 
GGA results in this Table are from non-selfconsistent calculations which employ the LDA charge density and geometry. 







LDA 




GGA 




LDA 


LDA 


11 


Geometry 


Es-A 


Ea-a 


Es-A 


Ea-a 


ds-A 


dA-A 


1 


On top of O 


1.46 




0.88 




1.89 




1 


On top of Mg 


0.45 




0.18 




2.51 




1 


Hollow 


1.01 




0.43 




2.07 




2 


O-Cu-Cu-0 (A) 


1.07 


0.93 


0.58 


0.83 


1.99 


2.25 


2 


O-Cu-Mg-Cu-0 (B) 


0.90 


0.77 


0.47 


0.72 


2.11 


2.34 


2 


Standing on O (S) 


0.96 


0.82 


0.71 


0.96 


1.87 


2.20 


3 


Line on O (A) 


0.96 


0.86 


0.52 


0.74 


2.02 


2.30 


3 


Triangle on O (B) 


0.85 


0.82 


0.45 


0.70 


2.06 


2.51 


4 


Rhombus on O (A) 


0.65 


1.08 


0.28 


0.93 


2.14 


2.32 


4 


Line on O (B) 


0.80 


0.89 


0.39 


0.81 


2.02 


2.59 



FIG. [l]. Free clusters obtained in LDA by relaxing their adsorbed geometries. Distances are in A . 
FIG. ^. Copper adsorption sites on MgO(lOO). 

FIG. ^. LDA potential energy surface for a Cu adatom on MgO(lOO). Energy units are in eV and taken with respect to the 
lowest lying state, i.e. the copper on top of the oxygen site. 

FIG. ^. Geometries of the supported Cu n clusters. For each case, we show both the starting geometry (top or left) and the 
optimized structure (bottom or right). 

FIG. [|. Dimer binding energy as a function of the angle 9 (see text for definition) as computed in the GGA. The line is only 
meant as a guide for the eye. 

FIG. ^j. Binding energies Ea-a and Es-a of copper clusters on MgO(lOO). Legends refer to the geometries of Fig. ^. GGA 
results are obtained through a non-selfconsistent method (see text for details). 

FIG. [?]. LDA density of states of free clusters. All DOS have been scaled down so that total charge is unity. 

FIG. ^. LDA density of states for the adsorbed Cu n clusters (n = 1 to n — 4 from top to bottom). All densities have been 
scaled down so that total charge is unity and aligned on the top of the oxygen 2s band. 

FIG. []. S, p, d projections of the DOS for a Cu adatom adsorbed at the oxygen site. 

FIG. 0. LDA charge densities for a Cu adatom on top of an oxygen site. This view is along the (100) plane, containing both 
oxygen and magnesium atoms of the slab. Crosses denote oxygen and copper atomic positions. 

FIG. [ll]. Density of states for a Cu adatom at the oxygen site calculated with different exchange-correlation functionals (bold 
line). All calculations are fully selfconsistent and include the optimization of the atomic structure. For comparison, the DOS 
of the clean surface (thin line) is included. The DOS have been aligned on the top of the corresponding clean surface's oxygen 
2s band. 
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